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The hindered diffusion model is introduced. It is a con- 
tinuum model giving the dynamics of a conserved density. 
Similar to the spin-facilitated models, the kinetics are hin- 
dered by a fluctuating diffusion coefficient that decreases as 
the local density approaches some geometrically constrained 
value where the mobility goes to zero. The model leads in a 
natural way to activated dynamics at low temperatures and 
high densities. In a well defined approximation the theory is 
shown to be compatible with mode-coupling theory. Indeed in 
the simplest form of the model we find an ergodic-nonergodic 
transition into a phase that supports activated kinetics and 
nonergodic behavior. 



We introduce a simple model for the structural glass 
transition that naturally leads to activated dynamics 
at low temperature and high densities. This model 
is continuum generalization of the spin-facilitated (SF) 
dynamics 1 ^ 4 that have been put forth as models for 
glassy dynamics. It is more physical than the SF mod- 
els since it is formulated in terms of a density field p 
that is conserved. As in the SF models, the equilibrium 
static behavior can be chosen to be noninteracting. In 
the continuum description this corresponds to choosing 
Gaussian static statistics. The model can also be used 
to make contact with mode-coupling theory 5 ' 6 . The key 
new ingredient in this approach is to choose a bare diffu- 
sion coefficient 7 that is density dependent D(p). Techni- 
cally this complicates the problem compared to the sim- 
plest field theoretical models that have been proposed 8,9 . 
This is because we must treat multiplicative noise. Gen- 
erally this has been a difficult proposition. However, 
there are certain constrained models, typically with con- 
servation laws, where the treatment 10,11 of multiplicative 
noise is tractable. The general class of models includes 
those discussed by Dean 12 , Kawasaki and Miyazima 13 , 
and Miyasaki and Reichmann 14 . The important differ- 
ence is that they deal with a bare diffusion coefficient 
that is linear in the density. We study here the situa- 
tion where there is a high density constraint on the bare 
diffusion coefficient. 

The model we study here is given in Langevin equation 
language by 



dp 
dt 



V 



V • (g P rf) 



where D(p) is the bare diffusion coefficient defined below, 
ff is Gaussian-white noise with variance 

(Va{x,t)r h (y,lf))=2kBTD Q 5 ce1 6(x-y)S(t-l/) , (1) 

and the effective Hamiltonian 



the multiplier g p — , j1 
Tip can be taken to be quadratic in p: 

U p = J d d xid d x 2 ^Sp(xi)x' 1 (x 1 -x 2 )Sp(x 2 ) 

where Sp(xi) = p{x\) — pa. 

The field theoretic version of this model corresponds 
to the MSR 15 action given by 



A = / <fxdt 



D(p)(yp) 2 + ip 



where p is the usual field conjugate to p. The theory 
from this point of view will be discussed elsewhere. 

We will analyze this problem using the Fokker-Planck 
(FP) formalism 16 . In this approach the density time cor- 
relation function is given by 

C(x 1 - x 2 ,t) = (<5p(x 2 ) e -^ t <5p(x 1 )) 

where the average is an equilibrium average over the 
static distribution e~ /3n, '/Z and the adjoint FP opera- 
tor is defined by 



D p = I d d x I d d y 



kuT- 



r P (x,y) 



5p(x) Sp(x) 
where 

r p (x,y)=V x -V y (D(p(x))S(x-y)) 
It is technically important 17 that 

d d z— T r p (z,x) = . 



s P (y) 



Sp(z 

We study a model with a simple physical choice for 



D(p) = D Q 9( Pc - p) 



(.Pc - p) 



Pc 



(2) 



where Do, introduced in Eq.(I), and p c are positive pa- 
rameters. This form reflects the geometrical fact that 



as the density increases it is more difficult for particles 
to move and there is some local density p c above which 
a particle is stuck. This is physically similar to spin- 
facilitated models where mobility is diminished if the 
local environment is blocked. We call this model the 
hindered diffusion model. More realistic models involve 
putting in more static structure and coupling the density 
to the other slow variables 18 in the system. 

In order to treat the high density constraint adequately 
we have organized the theory using the memory function 
formalism developed in Refs. 19 and 20. The Laplace- 
Fourier transform, 

C{k,z) = -iJ^ e"*y<f , a;ie ik -( Xl - Xa )C(xi-X2,t) , (3) 

satisfies the kinetic equation 



The activated dynamics results from the fact that if 
Po > Pc, then one can have motion only if there is a lo- 
cal fluctuation where p < p c . At low temperatures such 
fluctuations will be rare. 

Putting Eq.(5) back into Eq.(4), assuming 
and inverting the Laplace transform gives 



C(M) = e 



-k 2 D x -\k)t 



C{k) = e 



-r (fe)t 



C(k) 



(10) 



[z + K(k,z)]C(k,z) = C(k) 



(4) 



where C(k) — ksTx(k) is the static density correlation 
function and K(k,z) = K^(k) + K^(k,z) is the mem- 
ory function. We consider first the approximation where 
we keep only the static part of the memory function: 
K( s \ Formally is of order Do while the dynamic 
part K^ d \k,z) is of order D 2 . We will discuss approxi- 
mations for K^ d \ including mode coupling contributions, 
later. 

The static part of the memory function is given, with- 
out approximation, by 

K {s) {k) = il3- 1 k 2 DC-\k) =ik 2 D X ~ 1 (k) = iT (k) . (5) 

where D — (D(p)) is the average diffusion coefficient. If 
po = (p) is the average uniform density and S = ((Sp) 2 ), 
then, because the static fluctuations are Gaussian, we 
must have for the singlet probability distribution 



P[a] = (6(a-p(l)) 



(6) 



The parameter S is related to the static structure factor 

by 



S = 



d d k 



(27T) 



2k B T X {k) 



(7) 



and S is proportional to the temperature. The average 
diffusion coefficient can be evaluated as 



D = 



PqDq 



/ 7re(er/(e) + er/(e )) 



— e 



(8) 



where we have introduced the dimcnsionless parameters 



(Pc-Po) 

V2S 



1/Vr and e = -^t=. In the high density 

e < 0, low temperature limit |e| and eo are both large and 
we find that the physical diffusion coefficient is activated 
in temperature: 



D = D, 



Po 



1 



Or 2\pK£n 



1 

2^2 



4? + - 



(9) 



One has dramatic slowing down for e < and large. 

The dynamic part of the memory function is the sum 
of two pieces: 

KW{q, z)C(q) = (q; z) = F(q; z) + T sub (q; z). 

In coordinate space, the first piece is given by 

f (xi - x 2 ; z) = -(v( X2 )R(z)v( Xl )) (11) 

where the current is defined by w(xi) = iD p p(xi) and 
R{z) = [z + iDp]^ 1 is the resolvent operator. The Fourier 
transform of the second term, the subtraction part, is 
given by 

T sub (q:z) = (v^R(z)p (l )C- 1 (q;z)(p^R(z)v (l ) . (12) 

In carrying out the designated averages it is useful to 
write the current in the form 

v(l) = -iJ2^l f do-D{o)g„{l)p{l) . (13) 

where g a (l) = 6(a — p(l)) and we have introduced the 
field 

~p{\) =p ai (xi,ti) = W°l J d d x 2X - 1 (^i -x 2 )p(x 2 ,ti) 

In evaluating we make the simplest approxima- 
tion. We assume that the fields p appearing in the av- 
erages in Eqs.(ll) and (12), can be treated as Gaussian. 
Corrections to this approximation will be discussed else- 
where. We need, in the coordinate and time regime, the 
notation 

£7(12) - (6p(l)5p(2)) 



C{la u 2) = (p{l)5p{2)) 



C(lai,2a 2 ) - (p(l)p(2)> 

After doing the Gaussian averages, the subtraction term, 
Eq.(12), is given by 

r sub ( Xl -x 2 ;z) = Y1 V^V^C(x iai ,x 2 a 2 ;z)^ 2 . 



The more complicated contribution to the dynamic part 
of the memory function and is given, in the time regime, 



by T(12) = — (v(2)v(l)). After carrying out the average 
we find 



1 2 
x<5(q - ki - k 2 - k 3 )ki • q k 2 • q- (x~\k 2 ) - X^Oi)) 



f(i2)= V£V; 



+C(lai,2)C(2a 2 ,i; 



d 



dC(12) 



C(lai,2a 2 ) 



W(12) 



(14) 



where PF(12) is the diffusion-diffusion correlation func- 
tion that has the power-series representation 



W(12) = <£>(1)D(2)>=£ 



(C(12)) £ p 2 



(15) 



where we have introduced the average of the derivatives 
of the bare diffusion coefficient: 



d f 



(16) 



It is convenient to break VF(12) into three pieces 



Dr 



S ~ 



W{l2) = D 2 +D 2 C{^)+\ — ) ^-W(f) , (17) 

where the remaining sum given by W(f) can be rewritten 
as an integral 



f 

W(f) = / 
Jo 



sin 1 f 



dz(f — sin z)e <i+«™ *> 



(18) 



where /(12) = C(12)/S. 

We can associate distinct terms in with the three 
terms in VF(12) given in Eq.(17). The D 2 term can be 
shown to be equal to — f s „b(12) and cancels the subtrac- 
tion. The D\ term gives 



-.(d) 
MC 



^C(k,t)C(q-k,t) 



[q-[k X - 1 (fc) + (q-k)x- 1 (q-k)]] : 



(19) 



and we have the standard mode-coupling contribution. 
The term proportional to W in VF(12) gives rise to the 
contribution 



fa t) - a 2 I — 



n,\ 2 

2vr 



x X - 2 (fc)k-q/(fc,t)T^(q-k) 



where W(k) is the Fourier transform of W(f) given by 
Eq.(18) and W'(k) is the Fourier transform of -§jW(f). 
The dynamic part of the memory function is given by 
K^(q,t) = T ( ^(q,t)/C(q) and the static part by Eq.(5) 
These results go back into the kinetic equation given by 
Eq.(4). It turns out that we need to make some ad- 
ditional manipulations before looking for a solution for 
C(q,t). 

If we take the inverse Laplace transform of the kinetic 
Eq(4), we find 

C{t)-iK {s) C{t)- ( ds K {d \t-s)C(s)=0 . (20) 
Jo 

It is easy to see that Eq.(20) with K^ d \t) given by 
r< d )(t)/C*, does not lead to an ergodic-nonergodic tran- 
sition, instead numerical solutions blow up. Kawasaki 21 
suggested that the kinetic equation, Eq.(4), be rewritten 
in the form 



z + 



K (s) 



C{z) = C 



(21) 



1 + K(°)N(z) 

Comparing with Eq.(4) we can solve for N(z) to obtain 



N(z) 



(22) 



In a perturbation theory calculation where we treat 
K^ d \z) as small we have 



N(z) 



(23) 



So N(z) can be determined using perturbation theory. 
Next we note that Eq.(21) can be written in the form 

(1 + K^N(z))(zC(z) -C) + K^C(z) = (24) 

Taking the inverse Laplace transform gives 

C{t) = -T C{t) -T f ds N(t - s)C(s) (25) 
Jo 

where r = — iK^ sets the time scale. This equation 
produces a nonergodic phase for the standard mode cou- 
pling kernel. 

The last step in determining N(q,t) is to choose the 
static correlation function C(q). The simplest assump- 
tion is that the static correlation function is a constant 
up to a cutoff: xil) — \ f° r <Z < A and zero for larger 
wavenumbers. In this case we have 



'D^S 2 f d d h d d k 2 d d k 3 ~,,,. ffh T^{q,t) = ^-(rq 2 ) 2 S 2 ( ^f(k, t)/(q - k, 



t) 



+- 



M 2 (d V c2 



2tt 



S 2 J^0^k- q ,f(k,t)W(q-k,t) (26) 



where f(k,t) is the Fourier transform of f(r,t) = 
C(r,t)/S. r( d \q,t) is related to N(q,t) by Eq.(23), 
which takes the form: 

T^fo t) 

N(q, t) = = 9olMc(q, t) + 9iI B mc(q, *) (27) 

q^U z r z C 

where in three dimensions the couplings are given by 



2ttS 2 (D x 



C \D 



4tt5 3 (D 2 



C \D 



The functions lMC'(q,t) and Ibmc{<1^) follow from 
Eqs.(26) and (27) and the definitions of g and g\. Be- 
cause of the ratios of De the g's are not activated in 
temperature. 

The set of equations to be solved, Eqs.(25) and (27), is 
governed by the parameters p c , po, S and the wavenum- 
ber cutoff A. The static correlation function can be writ- 
ten as C(q) = 6tt 2 S/A 3 , which follows from Eq.(7). Time 
is measured in units of (Dr)^ 1 . 

In the dense frozen region where e < and is large 
m magnitude we have go = 3w and g\ = 3n . 
which suggests that for large e there will be an ergodic- 
nonergodic transition. 

We solved the coupled set of equations numerically. 
Let us first fix p c — 1.0 and the momentum cutoff A = 1.0 
and look at solutions C(q,t). If we begin with p$ = 0.5 
and S = 0.1 we find that C(q,t) decays exponentially 
with a decay rate, because of the conservation law, pro- 
portional to q 2 . Thus the kinetics are less sensitive to 
the cutoff as time evolves. If we increase S to 0.2 we 
see little change in the dynamics. However, it is easy to 
see that the couplings g n and g\ increase with S. While 
the model is only physically applicable 22 to structural 
glasses for small S, we can look at solutions for larger 
S and we find an crgodic-nonergodic transition (ENE) 
along the line Sh — Ah(ph — Po) 1 / 3 where Ah = 2.25 
and ph = 2.78. This holds even for po > p c . 

In the dense regime po > p c e < we find a line of 
ENE transitions for small values of S satisfying the rela- 
tionship S = A(p — p c ) 2 with a good fit for A = 0.245. 
Thus there is a regime where we simultaneously have 
activated dynamics and a nonergodic phase. As p n in- 
creases beyond p c , there is a value of po ~ Ph where the 
high S nonergodic branch meets the low S branch and 
for higher p n all states are nonergodic. One can look at 
the kinetics in the nonergodic regime and obtain a good 
fit to C{q,t)/C{q) = f(q) + A q /t. The ENE phase sepa- 
ration curves are determined by fitting to this form and 
choosing for fixed p and q = 0.4 that value of S that 
gives f(q = 0.4) = 0.5. 

Similar behavior is found for other choices of p c and 
wavenumbers. 
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